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1. INTRODUCTION 



EUGEN J. IONASCU 

Abstract. First, we calculate the Ehrhart polynomial associated to an arbitrary cube with integer 
coordinates for its vertices. Then, we use this result to derive relationships between the Ehrhart 
polynomials for regular lattice tetrahedrons and those for regular lattice octahedrons. These rela- 
tions allow one to reduce the calculation of these polynomials to only one coefficient. 

> 
O 

In the 1960's, Eugene Ehrhart (|13j.|14j) proved that given a d-dimensional compact simplicial 
complex in R n (1 < d < n), denoted here generically by V, whose vertices are in the lattice Z n , 
there exists a polynomial L(V,t) € Q[t] of degree d, associated with V, satisfying 

-i— > 

(1) L(V, t) = the cardinality of {tV} nZ", t€ 

It is known that 

> 

L(V, t) = Vol{V)t n + -VolidT)? 1 - 1 + ... + X {V) 

where x(7 ? ) is the Euler characteristic of V, and Vol(dP) is the surface area of V normalized with 
respect to the sublattice on each face of V . 

In [12], Eugene Ehrhart has classified the regular convex polyhedra in Z 3 . It turns out that only 
cubes, regular tetrahedrons and regular octahedrons can be embedded in the usual integer lattice. 
We arrived at the same result in [21] using a construction of these polyhedrons from equilateral 
triangles, and this is how we got interested in this line of questioning. This led us to the following 
simple description of all cubes in Z 3 . One takes an odd positive integer, say d, and a primitive 
solution of the Diophantine equation a 2 + b 2 + c 2 = 3d 2 (gcd(a, b, c) = 1). There are equilateral 
triangles in any plane having equation ax + by + cz = f, which can be parameterized in terms of 



two integers m and n. The side-lengths of such a triangle are equal to d\j2{ra 1 + mn + n 2 ). In 
order to rise in space from such a triangle to form a regular tetrahedron, one needs to have satisfied 
the necessary condition 
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Figure 1. Expand an equilateral triangle 



m 2 + mn + n 2 = k 2 for some odd k E N. 



Parenthetically speaking we remind the reader about Carroll Lewis ' conjecture on having infinitely 
many triples of Pythagorean triangles of equal area. The answer to that problem was to take 
solutions of ([2|) and construct triangles with sides 2uv, \u 2 — v 2 \ and u 2 + v 2 having (u,v) E 
{(m, k), (n, k), (m+n, k)}. This doesn't seem to be that strange of a coincidence, given the number 
of Diophantine equations that one need to solved in both cases. However, we could not find a 
reason of having ([2]) ciS ct sufficient condition. 

If ([2]) is satisfied, there are two possibilities. If k is a multiple of 3, then one can complete the 
triangle in both sides of the plane to a regular tetrahedron in Z 3 , and if k is not divisible by 3, 
then one can complete the triangle in exactly one side to form a regular tetrahedron in Z 3 (see 
Figure 1). Every such regular tetrahedron can then be completed to a cube in Z 3 with side- lengths 
equal to dk. Every regular octahedron in Z 3 is the dual of the double of a cube in Z 3 . We are going 
to make these constructions a little more specific in the last section. 

It is natural to ask the question that we think Ehrhart himself asked: "What is the form that the 
polynomial in (pQ) takes for these regular lattice polyhedra?" . The purpose of this paper is to answer 
this question for cubes (in a very simple way), and for regular tetrahedrons and octahedrons, the 
only such regular polyhedra in Z 3 . 

For completeness and due credit, we include Ehrhart 's idea in |15| to characterize all cubes in 
Z 3 . This was based on a theorem of Olinde Rodrigues: The set of 3-by-3 orthogonal matrices can 
be given by four real parameters a, b, c, d, not simultaneously zero, as follows: 



(3) 



a 



2 + b 2 + c 2 + d 2 



±1 



' a 2 + b 2 - c 2 - d 2 
2 (be - da) 
2(bd + ca) 
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It is clear that every cube in Z 3 can be translated so that a vertex becomes the origin and the 
three vectors defined by the three sides starting from the origin give an orthogonal basis for M 3 . 
Hence, one can construct a 3-by-3 orthogonal matrix from these vectors which has rational entries. 
Conversely, one can construct a cube in Z 3 from such an orthogonal matrix which has rational 
entries. In what follows we will do this association so that the vectors (points) are determined 
by the rows. The construction here is to take four integers a, 6, c and d in ([3]), simplify by 
whatever is possible and then get rid of the denominators to obtain the three vectors with integer 
coordinates that determine the cube. This construction is similar to the classical parametrization 
of the Heronian triangles. 

Our approach to the classification allows us to start in terms of the side lengths. However, 
Ehrhart's construction is useful in answering other questions about these objects. For instance, 
one can see that there are such cubes of any side length (other than the trivial ones, multiples of 
the unit cube) since every natural number can be written as a sum of four perfect squares. It turns 
out that there are only odd number side lengths for irreducible cubes, i.e. a cube which is not an 
integer multiple of a smaller cube in Z 3 . 

Let us begin with some of the smallest irreducible cubes. We introduced then here by orthogonal 
matrices with rational entries and defined up to the usual symmetries of the space (equivalent classes 
modulo the 48-order subgroup of all orthogonal matrices with entries or ±1, denoted by <S ) . As 
we mentioned before, and this will make a difference, the cubes are essentially determined by the 
rows. Obviously, the Ehrhart polynomials are identical for all cubes in the same equivalence class. 

We will denote the Ehrhart polynomial for an irreducible cube Ci of side- length £ = 2k — 1, 
k G N, by L(Ce,t). From the general theory we must have 

(4) L(Ce, t) = £ 3 t 3 + Ait 2 + A 2 t + 1, t e Z, 

where Ai is half the sum of the areas of the faces of the cube Ce, each face being normalized by 
the area of a fundamental domain of the sublattice contained in that face. The coefficient A2 is in 
general a problem (see, for example [5]) but in this case it takes a simple form as we will show in 
Section 3. 

For the unit cube C\ = I (the identity matrix), obviously, L(C\,t) = (t + l) 3 . There is only one 
cube (right or left equivalence classes modulo S a ) for each I = 2k — 1 for k = 1,2,3,4,5 and 6: 
Ci = I, 

-2 6 3 " 
3-2 6 , 
6 3-2 



C 3 := 



-1 2 
2 -1 
2 2 



4 


Co 







3 


-4 













5 
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Co 



7 4-4 
4 1 8 
-4 8 1 



, and Cii := — 



For k = 7 we have C13 := ^ 



-3 
4 
12 



12 4 
-3 12 
4 -3 



3 13 



2 
9 
6 



6 
-6 
7 



and an extra orthogonal matrix: 



5 
12 




12 
-5 
13 



One peculiar thing about the Ehrhart polynomials associated with these cubes so far, is that there 
is an unexpected factor in their factorization: 



L(C 3 ,t) = (3t + l)(9t 2 + l), L(C 5 ,t) = (5i + l)(25t 2 + 2t + l), L(C 7 ,t) = (7t + l)(49t 2 - At + 1), 
L(C 9 ,t) = (9t + 1)(81£ 2 - 6i + 1), L(Cn,t) = (lit + l)(121i 2 - 8t + l), 

L(C 13 ,t) = (13£ + 1)(169£ 2 - lOt + 1), and L(C 13 ,t) = (13* + l)(169t 2 + 2t + 1). 
This suggests that 



(5) 



L(C e ,t) = (£t + l)(£ 2 t 2 + at + l), t e Z, and some a £ 



We can easily prove that this is indeed the case for cubes of a special form, like C5 and C13 
above. Let us consider a primitive Pythagorean triple (a, b, c), with a 2 + 6 2 = c 2 . In the xy-plane, 
we construct the square with vertices 0(0,0,0), A(a, b, 0), _B(a — 6, a + 6,0), and C(— 6, a, 0). We 
then translate this face along the vector c k to form a cube of side-lengths equal to c. Let us denote 
this cube by C a ^ c . It is easy to argue that because we have a primitive Pythagorean triple (a, b, c), 
we have no lattice points on the sides of OABC, other than its vertices. The coefficient Ai in Q, 



is equal to 5 (c 2 + c 2 + 4c) because two of the faces have to be normalized by 1 and four of the faces 



i(c 2 + c 2 - 

have to be normalized by c(-) = c. By Pick's theorem, applied to OABC, we must have 

2 ^{points on the sides} . . r . . . 

c = \-#{mterior points of OABC \ — l = ^{interior points of OABC \ + l. 

Hence the number of lattice points in the interior of OABC is c 2 — 1. Therefore the number of 
lattice points in C a ^ c is (c+ l)(c 2 + 3) = c 3 + c 2 + 3c + 3. The polynomial L(C a ^ c , t) = c 3 t 3 + (c 2 + 
2c)t 2 + (c+2)t + l = (ct + l)(c 2 t 2 + 2t + l) satisfies exactly the condition L(C a>b>c , 1) = (c+l)(c 2 +3). 
We then have shown that ([5]) is true for infinitely many cubes Cg. 



a 
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Proposition 1.1. Given a primitive Pythagorean triple, a 2 + b 2 = c 2 , the cubes in the class of 
a b 

, have the same Ehrhart polynomial given by 



a,b,c 



-b a 
0c 



L(C C , t) = {ct + 1) ((ft 2 + 2t + 1) , t € N. 



The general formula is proved in Section 3. Section 2 is basically dealing with the second 
coefficient in (jl]). In Section 4, we look at the Ehrhart polynomial for regular tetrahedrons and 
regular octahedrons with lattice vertices. We show some nice relationship between the two and 
give a method of computing them in terms of Dedekind sums. 

2. The coefficient Ai 

Let us prove the following lemma which is more or less contained in |23| or it may already be a 
known result but we include it here for completeness. 

Lemma 2.1. Let a, b and c integers such that gcd(a,b,c) = 1. The sub-lattice C := {(x, y, z) G 
l?\ax + by + zc = 0} is generated by two vectors it and r such that \ u X r | = \/ a 2 + b 2 + c 2 . 

Proof: The statement is clear if two of the numbers a, b or c are zero. Otherwise the vectors 
it 



-(-6, a, 0), it 



-(— c, 0, a) and w 



g3(^)V - - g^d(^v — - - i^d(6^( >- c > 6 ) correspond to points in £. 

We first show that C is generated by u , it and ~w (a point P in V a< b tC is identified by its position 
vector OP as usual). Let (x,y,z) be an arbitrary point in £. If we define oo := gcd(a,b), then 
because gcd(a, b, c) = 1, we need to have z = loz' with z' 6 Z. Also, the existence of integers A: and 
/ such that ka + lb = u is insured by the fact uj = gcd(a, b). This means that we have 

(x, y, z) — z'[gcd(a, c)klt + gcd(b, c)£l$] = (a', ft , 0) G £. 

Since aa' + 6/3' = we see that (a',/3',0) = Alt for some A G Z 3 . This proves that the lattice L 
is generated by it, v and w. In fact, we showed even more: the the lattice C is generated by it 
and r := gcd(a, c)klt + gcd(b, c)€w, where ka + ib = uj. 

We observe that it x it = - gc ^ a c ) («~? + 6 J + c A; ) and similarly 17" x = - ge ^ fe ^ (a~/ + 
6 j + ck). Hence the area of the parallelogram determined by it and is equal to 

| it x 7*| = \ — (ak + bl)(a~t + b ~j + ct)\ = \J a 2 + b 2 + c 2 . ■ 

w 

Let us now assume that we have an arbitrary cube in Z 3 , 



(6) 



ai b\ c\ 
a2 62 C2 
03 b 2 c 3 
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with ai, bi and q integers such that cuaj + b$j + CjCj = 6ij£ 2 for all i, j in {1,2,3}. We define 
di := gcd(aj, 6j, Cj). It is clear that the di are divisors of i. Let us also introduce the numbers 
d\ = i = 1,2, 3. Then, we have the following expression for the first coefficient in @. 

Theorem 2.2. The coefficient Ai in is given by 
(7) A x = t{d x + d 2 + d 3 ). 

Proof: We use Lemma 12.11 for each of the faces of the cube. Opposite faces contribute the same 
way into the calculation. Say we take a face containing the points (ai,&i,ci) and (02,62,02). The 
irreducible normal vector to this face is clearly ^(03,63,03). The area of a fundamental domain 
here is given by ^ jp(a\ + 6| + c|) = d' 3 . By the general theory Ai = |(2y- + 2^- + 2y) = 
£(di + d 2 + d 3 ). ■ 
It is natural to ask at this point whether or not it is possible to have all of the dj's larger than 
one. It turns out that this is possible and as before, in our line of similar investigations, the first t 
with this property is I = 1105 = 5(13) (17): 

-65 156 1092 " 
420 1015 -120 . 
1020 -408 119 

Corollary 2.3. For a matrix Cg as in (EJ), such that Ci and C^ 1 are in the same equivalence 
class (left or right) modulo S 0) we have 



(8) 



C 



1105 



1105 



(9) d 1 +d 2 + d 3 = gcd(ai, a 2 ,a 3 ) + gcd(&i, b 2 , b 3 ) + gcd(c 1; c 2 , c 3 ). 

Proof. The Ehrhart polynomial must be the same for the corresponding cubes in the same 
equivalence class. ■ 
We believe that this corollary applies to all t < 1105, and of course to a lot of other cases, but we do 
not have a proof of this. A counterexample to the hypothesis of this corollary is given by the matrix 
given in ([8]). In this case, d\ + d 2 + d 3 = 35 and gcd(ai, a 2 , a 3 ) + gcd(&i, 6 2 , h) + gcd(ci, c 2 , c 3 ) = 7. 

3. The coefficient A 2 

The main idea in calculating the coefficient A 2 is to take advantage of the fact that every cube 
defined by (|6]) can be used to form a wandering set W for the space under integer translations 
along the the vectors Ti" = (ai,6i,ci), p = (a 2 ,6 2 ,c 2 ), and ~f = (03,63,03), i.e., 



o 

M 3 = (J (W + tct + j^ + k^), 
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Figure 2. Wandering set determined by the cube 



where (J means a union of mutually disjoint sets. 

The wandering set W we will consider here, associated with a generic cube as in Figure [21 is the 
set of all points formed by the interior points of the cube to which we add the points of the faces 
OAEC, OADB and OBFC except the (closed) edges AD, DB, BF, FC, CE, and EA. It is easy to 
see that such a set is indeed a wandering set. In our setting we think of of, p and as the vectors 
OA, oii and OC. With the notation from the previous section, we have the following result. 

Theorem 3.1. The coefficient A 2 in §4ty is equal to d\ + di + c?3 • 

Proof. Let us denote by k the number of lattice points in W. For n £ N, then the number of 
lattice points in 

o 

\J (W + + k^), 

M,fce{l,2,...,n} 

is equal to rfik. On the other hand this number is equal to L(Ci,n) + K, where K is the 
number of lattice points on three big faces of nC. It is easy to see that K is 0(n 2 ) and so 
k = \{L{Ct,n) + 0(n 2 )) = £ 3 . We now need the following result which is well known. 

Theorem 3.2. (Tlhrhart-Macdonald Reciprocity Lawj Given a compact simplicial lattice 
complex P (as before) of dimension n, then 

L(P,t) = (-l) n L(P,-t), t e N. 

Hence, according to this theorem the number of lattice points in the interior of Ci is £ 3 — Ai + A2 — 1. 
So, the number of lattice points on the boundary of Q is 2Ai + 2. Let us denote by a the number 
of lattice points in the interiors of the sides OA, OB and OC. Then we must have 



2Ai + 2 = 2[k - (£ 3 - Ai + A 2 - 1)] + 2a + 6 =>■ A 2 = a + 3. 
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Since, a = (dx — 1) + (d 2 — 1) + (d 3 — 1) the claim follows. ■ 
Putting these facts together we obtain. 

Theorem 3.3. Given a cube Ci constructed from a matrix as in {6j), its Ehrhart polynomial 
is given by 

(10) L(C e ,t) = {£t + l)[ft 2 + {d 1 +d 2 + d 3 -£)t + l], t e Z. 

There are some natural questions at this point. One of them is: "What is the maximum number 
of lattice points that one can contain in a lattice cube of side lengths £?" We have the following 
corollary to the above theorem. 

Corollary 3.4. Given a cube Cg constructed from a matrix as in ^) , the maximum of lattice 
points inside or on the boundary of this cube cannot be more than (£ + l) 3 . This value is attained 
for the cube tC\. 

Proof. Since di is a divisor of £ we must have dj < I and so the corollary follows from (110p . ■ 
What is the answer to above question, if we restrict the problem to irreducible cubes? It is clearly 
a more complicated problem and it depends heavily on I. 

Finally, another question which seems natural here is: "What is the cardinality of the set of all 
Ehrhart polynomials associated to irreducible cubes of sides £1" 

4. Regular Tetrahedrons and regular Octahedrons 

We remind the reader that the cube in space (Figure [3]) is determined by an orthogonal matrix 
as in ([6]) by taking its vertices O (the origin), A, B, C, D, E, F and G, whose position vectors are 
ol = -ct= ( ai ,6 1 ,c 1 ),0^ = ^ = (a 2 ,b 2 ,c 2 ),od = t = (a 3 ,b 3 ,c 3 ),OTJ = ^ + ^,Op = ^+^, 
0% = -f + tt, and 0& = tt + ~$ + -f. 

In [21| . we rediscovered Ehrhart's characterization ([E]) of all regular polyhedra which can be 
imbedded in Z 3 . Only cubes, tetrahedrons and octahedrons exist in Z 3 and there are infinitely 
many in each class. We have constructed all of these out of equilateral triangles. In general, once 
a tetrahedron is constructed, it can always be completed to a cube. Vice versa, for a cube given 
by ([6]) , there are two regular tetrahedrons which are shown in Figure O They are in the same 
equivalence class modulo the orthogonal matrices with entries ±1, denoted earlier by Sq. Regular 
octahedrons can be obtained by doubling the coordinates of the cube Ci and then taking the centers 
of each face. This procedure is exhaustive. An octahedron in the same class can be obtained by 
simply taking the vertices whose position vectors are ±"3", ± and We will going to use the 
notations Ti and Oe for the tetrahedrons and octahedrons constructed in this way from Ci. If we 
deal with irreducible such objects, i.e. they cannot be divided by an integer to obtain a strictly 
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Figure 3. A cube determines two tetrahedrons 



smaller similar polyhedron, then we may assume that I is odd. The Ti and Oi have side-lengths 
equal to Ey2. From the general Ehrhart theory (see [3]) we must have 



11 



L(Tt, t) = —t 3 + fiit 2 + fi2t + 1, L(O e ,t) 



U 3 



t 3 + v x t 2 + u 2 t + 1. 



Let us first look at some of the examples with the smallest side- lengths. 

with 





' 1 


1 


" 




' 1 





" 


Ti : = 





1 


1 


, O l := 





1 







1 





1 










1 



For the next sides, 



*) 


t 3 




U 






= 


4 „ 


„ 8t 


3 


- + 


t 2 + 


y 


+ 1, and L(Q\, 


-f + It H 

3 3 






" 1 


1 


4 " 




' -1 


2 


2 " 




T 3 




1 


4 


1 


, O s := 


2 


-1 


2 


, with 






4 


1 


1 




2 


2 


-1 





L(T 3 , t) = 9t 3 + -t 2 + + 1, and L(0 3 , t) = 36t 3 + 9t 2 - t + 1. 
These polynomials were computed initially with the help of a computer. 



4.1. Coefficients [i\ and v\. From the general theory we know that these coefficients can be 
computed in terms of the areas of their faces normalized by the area of the fundamental domains 
of the sub-lattice of Z 3 corresponding to that face. Since the number of faces for Ot is twice as big 
as for Ti and since the parallel faces are in the same equivalence class, we must have v\ = 
Let us introduce the divisors 
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(12) 



Dx = gcd(ai + a 2 + a 3 , h + b 2 + 63, c x + c 2 + c 3 ), 
£> 2 = gcd(ai + a 2 - a 3 , &i + 6 2 - 03, ci + c 2 - c 3 ), 
Z>3 = gcd(ai - a 2 + a 3 , 61 - 6 2 + 63, ci - c 2 + c 3 ), and 



£>4 = gcd(-ai + a 2 + a 3 , -61 + o 2 + 63, -c\ + c 2 + c 3 ). 

Observe that the vectors ~ct + + , ~cl + ~$ — ^f, ~ct — p — 7jf + ~$ + are vectors normal to 
the faces of the Tg. By Lemma [2. 11 we see that the area of each fundamental domain corresponding 
to a face of Tg is given by one of the numbers 

Proposition 4.1. The coefficients [i\ and v\ in Ul\) . are given by 



(13) 



Ml = "J = «Pi +^2 + ^3 + £> 4 )/4. 





" 7 


-1 


" 




' 4 


CO 


" 


T 5 : = 


4 


3 


5 


, 5 := 


3 


-4 







3 


-4 


5 










5 



This explains the next examples which were obtained by brute force computations with Maple: 



, with 



L(T 5 , t) = ^t 3 + 5t 2 + \t + l, and L(0 5 ,t) = ^ 3 + Wt 2 + ^-t + l. 

4.2. Coefficients u 2 and v 2 - Let us observe that the cube in Figure [3] can be decomposed into 
four triangular pyramids OABC, DABG, FCGB, and EGCA, which can be translated and some 
reflected through the origin to form half of Oi, and the regular tetrahedron Tg. We remind the 
reader of some notation we used in the proof of Theorem 13.11 we denoted by a the number of 
lattice points on the interiors of the edges OA, OB, and OC. We showed that a = d\ + d 2 + 03 — 3. 

Let us balance the number M of interior lattice points of C( using the above decomposition. 
According to Theorem 13.31 and Theorem 13.21 we have 

) = (£ - i)(e 2 - - d 2 - d 3 + e + 1) = e - (di + d 2 + d 3 )(£ - 1) - 1 



M 



-L(C e 



Some of the lattice points counted in M are in the regular tetrahedron, and are counted by 
L{Tt, 1) = y + Hi + n% + 1. From these we need to subtract the number of lattice points on 
the interiors of its sides, which we will denote by r, and also subtract 4 for its vertices. The rest 
of the points counted in M, are in the interiors of the four pyramids. If we multiply this number 
by two and add the number of lattice points in the interiors of the cube faces less r, we get the 
number of interior points of Ot minus 2a + 1. The number of interior lattice points of the cube 
faces is equal to 2Ai + 2 — 4cr — 8. In other words, we have 
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£ 3 4/ 3 
2(M - — - m - n 2 - 1 + r + 4) + 2Ai + 2 - 4a - 8 - t = — - v x + v 2 - 1 - (2a + 1). 

Taking into account that z^i = 2/xi and Ai = (di + d 2 + d^)£ = (a + 3)£, this can be simplified to 



(14) i/ 2 + 2/x 2 = 6 + r. 

We close this section concluding what we have shown. 

Theorem 4.2. For a regular tetrahedron Tj? and a regular octahedron Op constructed as before 
from an orthogonal matrix with rational coefficients as in ([6}), the coefficients \x 2 and v 2 in 
satisfy 



(15) u 2 + 2fi 2 = (d 1 +d 2 + d 3 + d 4 + d 5 + d e ), 

where d\, d 2 , d% are defined as before and d 4 = gcd(ai — a 2 , b\ — b 2 , c\ — c 2 ), d§ = gcd(ai — 03, b\ — 
h,ci - c 3 ) ; and d 6 = gcd(a 3 - a 2 ,6 3 - b 2 ,c 3 - c 2 ). 

We have tried to find another relation that will help us find the two coefficients but it seems like 
there is not an easy way to avoid, what are called in [3], the building blocks of the lattice-point 
enumeration, the Dedekind sums. This numbers require a little more computational power and we 
are wondering if a shortcut doesn't really exist. One would expect that the answer to our questions, 
for such regular objects, is encoded in the coordinates of their vertices in a relatively simple way. 
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